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Introduction 


The primary goal of this project is to understand the sources of the near-Earth ambient magnetic field 
as observed by recent spacecraft surveys and surface variational magnetic observations so as to determine 
(as well as the data allow) the electrical properties of the crust and upper mantle. Also included in the 
study is the structure and changes on a short time scale of the core field which must be separated and 
identified in order to achieve the pirmary goal. 

The Magsat data collection interval (Nov 1979- April 1980) provides an opportunity to compare the 
vector field projections of ionospheric currents computed from surface data above the ionosphere as does 
the POGO data (1965-1971) for scalar projections. The limitation of Magsat is its sun-synchronous orbit, 
which only sampled low latitudes at dawn and dusk, whereas POGO, though only making observations 
of the scalar field, sampled all local times. Magsat operated at a lower altitude than POGO (down to 350 
km) whereas the orbits of the three POGO spacecraft ranged up to 1500 km and were never lower than 
about 400 km.— 


The planned program of this study is proceeding in parallel along several lines: 

1. Obtain the spherical harmonic coefficients needed to interpolate observatory hourly values. This is 
initially being done using observatory data taken during the 12 month interval July 1, 1979 through 
June 30, 1980 to span the Magsat data collection interval. Earlier data, taken during the POGO data 
collection (1965-1971) have also been used and will be part of the analysis following the Magsat 
work. The intent is to produce a set of hourly coefficients for the interval from about 1960 through 
1980 for these and future studies. 

2. Obtain an accurate high degree model of the static field so as to be able to compare the projected 
results from (1) with the deviations from Magsat to investigate how these potential expressions portray 
the field variations above the ionosphere. (We already know of deviations due to meridional and field 
aligned currents and plan to add the appropriate non-potential functions to allow their modelling.) 
Study and understand the ionospheric/magnetospheric interactions as are appropriate to observatory, 
Magsat and POGO data sets. 

3. Develop functions that will relate the external variations imposed on the Earth to its electrical 
response. 

4. Invert the results of the intemal/extemal relations to determine physically meaningful conductivity 
models of the crust and upper mantle. 

5. Study the short period secular variation determined from the interval 1960-1980 in regard to such 
topics as frozen flux and energy transfer as may be pertinent to a fuller understanding of the physics 
of the core dynamo. This requires separation of apparent core secular variation due to long term 
changes of externally imposed flux from that arising in the core. 

6. Look for the physical bases for the crustal signals for the range n = 14 - 50; whether they can be 
explained completely in terms of crustal magnetization. Our speculation is that there is the possibility 
that potentials driving weak currents through conductivity channels in the crust and upper mantle may 
be contributing to the crustal signatures. 
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Data 


The observatory data being used is indicated in Table 1. These are the data for which a complete year 
are available (July 1, 1979- June 30, 1980), and in hand. In addition to this set, there are others for which 
a partial year is available or more that has been found recently but not yet acquired or processed. These 
are noted in table 2. 

There are a number of stations for which there should be data, because they were collected by this 
project or entered onto the data base from printed listings. The reductions were only being done for data 
from September, 1979 through June of 1980. These had been turned over to the data center in 1985 or 
1986. It is thus unclear and subject to further investigation as to what happened to the data sent to them 
at that time. Figure 1 shows the locations of these stations. 
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Table 1: Magnetic observatories for which a full year’s data (July 1, 1979- June 30, 1980) are in hand 
and either adjusted if discontinuities, or otherwise processed into an analytic format (Most of these can 
be identified by reference to several listings provided by the World Data Centers for Geomagnetism or 
in such references as Langel (1987)). 


AIA 

GWC 

MNK 

AMT 

HER 

NUR 

AQU 

HON 

NEW 

BNG 

IRT 

NVS 

BOU 

KAK 

PMG 

BRW 

KGL 

SGE 

CBB 

LOV 

SIT 

CMO 

LRV 

SJG 

CNB 

MBC 

STJ 

CZT 

MEA 

TUC 

FCC 

MMB 

VIC 

GUA 


WIT 



YKC 


Table 2 Data which are available from sources other than the World Data centers or which are received 
but not yet validated or entered into the analytic base. Those in parentheses are available from other 
sources and are reported to have been validated, those in square brackets are not yet received and are in 
some raw form. The numbers following the brackets indicate the year and the number of element months 
reported to be available for that year. 


[ALE 79: 30] 

BLC 79 
BJI 80 

BLC 80 [BLC 79:33] 
CAN 80 [CAN 79:12] 
CNH 80 

GNA 79 [GNA 80:18] 
DRV 80 [DRV 79:43] 
[ebr 79:2] 

ESK 80 

FRD 80 (FRD 79) 

[FUR 79:11, 80:10] 

GNA 80 
GZH 80 
HAD 80 (HAD 79) 
HRB 80 

HUA 80 [HUA 79:12] 

ISK 80 

JAI 80 

KGD 80 

(KNY 80) KNY [79:36] 

(LER) 

MBO 80 [MBO 79:12] 
MIR 80 
MIZ 80 


(MOS 79) 

MMK 

80 [MMK 79:3 

LMM 

80 

LNP 

80 

LZH 

80 

(NGK 80) 

NKK 

80 

ODE 

80 

OTT 

80 [OTT 79:8] 

POD 

80 

PPT 

80 [PPT 79:12] 

RES 

80 [RES 79:8] 

RSV 

80 

SNA 

80 

SOD 

80 

SSH 

80 

SVD 

80 [SVD 79:35] 

TKT 

80 

TOL 

80 [TOL 79:12] 

[TOO 79:18] 

TRD 

80 

(VLA 79) [VLA 80:36] 

VSS 

80 [VSS 79:12] 

WHN 

80 

WHS 

80 [WHS 79:4] 

(YAK 79) 
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Theory 


The classical method of representing diurnal variations of observatory hourly data has been to compute 
an equivalent sheet current system at about 100 km altitude expanding the potential function of the 
observed variations in low degree spherical harmonics ( Chapman , 1919). This technique has been applied 
mainly to averages taken over quiet intervals after removing a base value for each station taken. Such 
removal implies that the equivalent electrical currents obtained are those relative to that at some local time 
(e.g. midnight), when such currents should be weak or non-existent. The usual coordinate system for 
making such computations has been latitude (or magnetic latitude) and some form of local time. A more 
modem look at such currents was by Suguira and Hagan (1967) who were able to compute the potential 
of the variations for each hour of Universal time for two selected 5-day intervals. 

The potential function of the time variational field was expressed by: 

4 n 

y(r,a f A,T) = a ^ cosmA + An ^ CT3 sin mAj + 

r»=] m= 1 

cos mA + sinm*] (j]" +1 ] PfCcotSl 

The variational field B is obtained from such functions by the negative gradient B = - W. 

A problem with determining the coefficients is the large gaps in spatial coverage of the data, even 
limiting the degree to n = 4. Past experience with such analyses has been that the least squares process 
would sometimes not converge, or would converge on very spurious results with too sparse a data 
distribution. It has often been found useful to supplement the distribution allowing the data for a given 
hour to enter the computation several times. One method used 3-hour running averages by considering 
time differences to be equivalent to longitude differences within each 3- hour interval. If the computation 
was done to produce values at time T, a value of the field difference for time T-l and T+l was also 
entered with longitudes -15° and +15° respectively from the longitude at time T. This is equivalent to 
using a box car filter with 3 point overlaps in the time domain. 

Our approach here is to develop the coefficients per UT day by expansion of the time terms in the 
coefficients. However, difficulties in obtaining detailed results, or, in some instances, results that will 
converge at all, may force us to adopt similar expedients as above. Also, because data will be obtained 
above the presumed source of ionospheric currents as well as below, we separate out the different 
potentials as noted in our earlier report as follows: 


V = Vi + v e s + Vi + v e m 

where the potentials are identified with the following sources: 


yj : core field 
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Vg : E/F layer ionospheric currents 

V' f : Internal Earth currents induced by external flux changes 

Vm : Magnetospheric currents external to the spacecraft shell as well as the Earth 

For our purposes here we expand the coefficients of the internal field into linear time derivatives so 
as to take into account secular variation only over a short period. This expression is thus: 


\ n+j * 

V'c — a ^[f] y {[ 9n + 9n C SO ] cos mcj> + [ h ™ + h™ C SO ] sin rruf) P™ CflD 


n= 1 


m=0 


The present limit to n c is 50 for the spatial terms and 7 for the first time derivatives (SV) as 
determined for the Magsat interval. This is the most recent model determined using least squares only for 
the high degree terms and will be published soon ( Cain et al. 1990). The secular change terms are those 
derived by Cain et al. (1989a) using a combination of observatory secular change and Magsat. This SV 
is not adequate and will need to be revised after the external field computation is complete because it is 
likely that some of the seasonal changes from ionospheric currents and their induced analogs are 
incorporated in the first time derivatives. 


The form of the Vfn terms is: 


Kn = T ^<7j e cos 0 + [j7j e cos <j> + /i] e sin <f>~\ sin 9 } 


The potential for the currents induced in the Earth from the changes in all external currents, both 
ionospheric and magnetospheric, is: 


n n 

V 1 , = a ^ f~"| ^ [G™ cos m<j> + H%' sin m9~] 

n— 1 m=0 

and for the ionospheric currents themselves: 


rt , n 

V s , = a y y [ G™cos rruf) + H% e sin ro<ft] 


n= 1 m=0 
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where the terms for G™ e and H™ e are: 

n 

gt = £(«£“«;« * ^/siny/) 

;= i 


;=i 


with expressions for the internal functions corresponding to these by replacing "e" by "i." 

The above are for computations below the ionosphere (i.e. at the Earth’s surface). However, above 
the ionosphere, the potential for the ionospheric currents changes so that, for a continuous vertical field 
across the sheet, the potential becomes an internal one: 


« = -ns 


C2n+l] f 


n 


i] [?] c ” i} xj 


sin 


mfi] 


It is clear that there is no way from surface data to distinguish the n = 1 tenms from the ring current 
from an n = 1 term from the ionosphere. Thus for the determinations prior to addition of satellite data, 
the adjustment for the potential Vfn is being turned off and left entirely to the potential V$- (Most other 
workers have ignored this term). 


For our studies we will make no selection of data for quiet, but instead compute the potential for all 
UT hours under consideration. Although the field aligned currents and the near equatorial meridional 
currents ( Yanagasawa and Kono, 1985; Cain et. al., 1989) are not yet added to the fitting functions, most 
of these effects are non-potential and are ignored initially. A reason why it is valid initially to ignore these 
currents is that the vertical component are non-potential and thus do not affect the potential function being 
developed. Their horizontal components will be potential, and will vary with altitude and affect the Magsat 
data especially. In polar regions the field aligned currents are near vertical and most of the main field 
analysis has used data when they were very weak. Recent work by Langel, Rajaram, and Purucker (1990) 
on the structure of the meridional currents near the dusk equator will be helpful in further modeling for 
these currents. 
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Potential field testing 

As a means of testing this program, some of the previously reduced data during the POGO interval 
were used. Data were from 34 magnetic observatories for the period Sept 22 - 27, 1967. By using the 
expedient of computing the variational potentials from the rediduals of the data from the average night- 
time values (21-2 hours local time), the coefficients were obtained as seen in Table 3. These appear to 
have about the correct amplitudes compared with those of Sugiura and Hagan (1967) (after converting 
theirs to Schmidt normalization). This can be especially seen in the large diurnal term for the n = 2, m 
= 1 external coefficents. 
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Table 3a 
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Table 3b 
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Another way to view such results is to look at maps of either the current function or the potential 
function. The two functions are similar. As noted earlier, the potential above and below the assumed 
ionospheric current sheet are different and can be also expressed in terms of the potential of a surface 
current as follows: 


V — p~] S n , for r)R or 

n = 1 

CO 

R ^ Cn+lD 

n=l 


'*■ 

r 

Jt. 


Sn, for r(R 


where S is: 


n 

S n = ^ cos m<f> + rjn sin m<p') 


m=0 


The current density J of such an ionospheric sheet can be shown (within an arbitary constant) to have 
components in the south and east directions respectively of: 


j, = -A 

9 2x 

CO 

y"f2n + r 

| lain 0^ 

1 

as„ 

df> 

ii 

CO 

C2 n + ID 

dS 7 

d$ 


The potentials for V,‘ can then be matched up to those for S H with the result that the contours of J and 
of are parallel. We here plot the potentials of V t = V‘ + Vj in order to visualize whether the program 
is producing reasonable patterns to also correspond in structure to those of equivalent ionospheric current 
systems. That is, the current direction is indicated by that of the contours and the current density by the 
crowding of the contours. 

Figure 2 is a representation of the contours of the sum of the potentials V t for 4 different universal 
times. There is clearly some problem with a discontinuity on the prime meridian due likely to the method 
of analysis. Otherwise, the cellular structure is not dissimilar from those plotted previously. In general the 
strongest currents are near the noon meridian with the northern focus earlier in time than the southern. 

However, there are sufficient differences from prior current patterns to make it imperative that these 
computations be scrutinized in detail and further tests made using synthetically produced data to assure 
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that all programs are producing accurate results. Also, instead of attempting to validate the POGO data 
tfiat went into this computation, which would require re-acquiring surface data for this interval, we plan 
instead to first concentrate on the Magsat interval. 

Internal/External Relations 

The development of a connecting matrix to express the response of the Earth to an imposed external 
field was developed by a prior NASA contract (NAS 5-9319) but never implemented completely due to 
the complexity of the functions required and the limitations of computing equipment at that time. The 
development that follows is thus an outline of the technique that might be a useful tool for analysis 

simultaneously to providing the basic functions needed to determine the Earth’s internal conductivity 
structure. } 


The premise was that there exists an impedance function z(9,<J>,<o) that relates the tangential electric 
field E and surface current density by the equation E = z J lan . That is, the solid Earth would be represented 
by a thin shell. The usefulness of such a function becomes apparent especially when one looks at the 
power spectrum of the diurnal variations (e.g. Parkinson, 1983, P. 261). These show that the major power 
is in harmonics and subharmonics of the diurnal period. Also, luni-solar harmonics are imposed from the 
lunar tidal terms interacting with the changing electron density. 

Thus it may be practicable to determine such a shell function at specific frequencies which represent 
most of the power in the quiet diurnal field patterns. There is some evidence that the disturbance field may 
also be so broken up into a small number of discrete frequencies. However, the theory does not require 
that a small number be used, and there may be ways of obtaining band passed averages if the spectrum 
does not have sharp peaks as it does during quiet conditions. 


It can be shown that if 


a - ( 2n+1 ) S in „ 

n = 1 
oo 

b = -R/ n Cn+13 C2n+l3 S- 

/ i in n 

n = l 

oo 

c = *£ (A*. + 4J 

n = l 

equations can be set up of the type 


Va*Vz + bz = c 

which should yield two real functions for each frequency given the proper variational input to be able to 
produce current functions of the type S n as expressed earlier. 
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Crustal Field 


The high degree field model from Magsat and surface observatory secular change data was initially 
developed by a series of papers now in print ( Cain et al. 1989a,b). This model was recently tested under 
this project and found to have questionable error distributions so was recalculated with adjustments to n 
= 50. This paper is to be published in September and a copy is attached herein. Also included in this copy 
are some added error distributions that did not appear in the printer version of the paper. 

It is not completely certain why the integral method applied to the Z component for obtaining the 
coefficients to n = 63 gave such poor results. Our supposition at this time is that it was due to the fact 
that the field above the ionosphere was seriously affected by the external currents. It was noted in the 
development of the model for correction of the "ionospheric currents" ( Cain 1989a) that the residuals were 
not uniform as to epoch and altitude. There were systematic biases depending on how the data were 
selected. However, due to the need to obtain results with less than adequate resources, we were unable 
to pursue these questions in detail. It would seem that the new model is adequate to use as a base for 
comparing the results of the surface model projected ionospheric current changes to satellite altitude. 

Future plans 

In the last progress report it was noted that we were faced with considerable uncertainty as our 
supercomputer company (ETA) had suddenly gone out of business. At this time we have now converted 
to an operational Cray YMP and also have an added Connection Machine 2. Our Sun workstation has also 
had its memory expanded to 24 mbytes and disk expanded to 900 mbytes. Thus such problems are now 
essentially behind us. 

Validating and completing the observatory hourly value set for the Magsat interval is high on the 
priority list. This will require more detailed looking at the baselines and comparing with averages 
computed by Cain and Kluth (1984, 1985). 


The procedure of computing the external and internal field variations from surface and satellite data 
needs to be rechecked using synthetic data. 

Further work needs to be done on the theory of relating internal and external source functions. 
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Appendix A - Numerical Experiments in Geomagnetic field modelling 
J. C. Cain, B. Holter, and D. Sandee 

Journal of Geomagnetism and Geoelectricity, 42, 961-972, 1990. 


This is the version as submitted to the above Journal. Several small errors were corrected in the final 
version which is not yet available. 
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